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ABSTRACT 

An understanding of the oscillations of differentially rotating systems is key to many 
areas of astrophysics. It is of particular relevance to the emission of gravitational waves 
from oscillating neutron stars, which are expected to possess significant differential ro- 
tation immediately after birth or binary merger. In a previous paper we analysed the 
normal modes of a simple system exhibiting differential rotation. In this complemen- 
tary paper we address the initial value problem for the same simple model using both 
analytical methods and numerical time evolutions. We derive a necessary and suffi- 
cient condition for dynamical shear instability. We discuss the dynamical behaviour of 
the continuous spectrum in response to an initial perturbation, and show that certain 
singular solutions within the continuous spectrum appear physically indistinguishable 
from the discrete modes outside the continuous spectrum. 
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1 INTRODUCTION 

The oscillations of differentially rotating systems are im- 
■ portant in many areas of astrophysics. They are of par- 
ticular relevance to neutron star astrophysicists, because 
, non-axisymmetric oscillations of compact stars are promis- 
ing sources of gravitational waves. Detection and interpre- 
tation of gravitational waves from oscillating neutron stars 
could answer many questions about neutron star physics 
that are difficult to answer using electromagnetic obser- 
vations, for example regarding the supranuclear equation 
of state. Emission of gravitational waves from neutron 
star oscillations will be strongest if the oscillations can 
grow via some dynamical or secular instability. The most 
promising instability scenarios arise for newly-born, accret- 
ing or newly-merged neutron stars. But it is in precisely 
these situations that we expect the stars to be differen- 
tially rotating. The latest simulations of core collapse in- 
dicate that the nascent neutron star emerges rotating differ- 
entially (Diminc lmeic r. Font fc M iillor 2002; Akivama et al 
I2OO3I: lOtt et ail [2to3). iFuiimotol J 19931) has shown that 
rapid accretion from a companion could drive differen- 
tial rotation in the surface layers of the neutron star, 
and a similar effect may occ ur during rapid accretion 
of su pernova fal l back material iWatts fc Andcrsson 2002). 
IShibata ~fc Urvul (|2p0(l) and lBaumgarte. Shapiro fc Shibatal 
l|20000 have shown that a long-lived massive differentially 
rotating neutron star may be generated as the result of 
the merger of two neutron stars. In addition, recent work 



on r-modes in neutron stars suggests that the modes may 
drive a uniformly rotating star into differential rotation via 
non-linear effects fRezzoUa et al"200rt 'Levin fc Ushomirskvl 
2001; Lindblom, Tohlinc & Vallisncri 2001). Understanding 
the effect of differential rotation will be critical if we are to 
model gravitational wave emission from these systems cor- 
rectly. 

Differential rotation introduces two interesting phe- 
nomena. F]irstl^^j_J^jna^MfiadtodTO insta- 
bilities (iPapaloizou fc Pringld Il984l . Il985l : iBalbinskil Il985l : 
lLuvten|[l990ir Secondly, the dynamical equations of the in- 
viscid problem are complicated by the presence of corotation 
points. These are points where the pattern speed of a mode 
matches the local angular velocity, and at these points the 
governing equations are formally singular. The range of fre- 
quencies for which modes possess corotation points is termed 
the corotation band. It gives rise to a continuous spectrum, 
the dynamical role of which is unclear. Neither issue has re- 
ceived sig nificant attention i n the stellar perturbation liter- 
ature fsee lWatts et all i2003h . hereafter referred to as Paper 
I, for a summary of previous work). 

We have begun to address the problem by modelling a 
simple system, a differentially rotating incompressible spher- 
ical shell. This model is simplistic, but exhibits the main 
features of the more complex three-dimensional problem. In 
Paper I we conducted a normal mode analysis for the shell. 
In the uniform rotation limit the shell problem admits only 
the standard r-modes as solutions. As the degree of differen- 
tial rotation increases, the frequencies of these modes grad- 
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ually approach the lower edge of the expanding corotation 
band. Most modes approach the boundary tangentially and 
appear to cease to exist when the degree of differential ro- 
tation exceeds a certain threshold value. In only a few cases 
do modes cross the corotation band boundary. 

For real frequencies in the corotation band the shell 
problem admits a continuous spectrum of solutions with dis- 
continuous first derivatives. In general, the continuous spec- 
trum solutions have both a logarithmic singularity and a 
finite step in the first derivative at the corotation point. 
For certain special frequencies, however, the step in the first 
derivative vanishes. We refer to these solutions as zero-step 
solutions. Since these solutions appear when regular normal 
modes cross into the corotation band, they are likely to be 
of particular importance. Furthermore, in difi'erentially ro- 
tating cylinders, these solutions can merge to give rise to 
dynamical instabilities (Balbinski 1985). Even though we 
do not observe such behaviour in our problem it is an in- 
teresting possibility that should not be overlooked in future 
studies. 

The character of the modes that cross into the corota- 
tion band varies. Only in one exceptional case, for the sim- 
plest rotation law examined, do we observe a mode crossing 
the corotation band boundary and continuing to exist as 
a regular stable mode. We think that such regular modes 
within the corotation band are unlikely to exist for more re- 
alistic rotation laws. In general, at all points where modes 
cross the corotation band boundary we observe the develop- 
ment of zero-step solutions. At some points we also observe 
the development of dynamically unstable modes. A neces- 
sary (but not sufficient) condition for dynamical instability 
is that the equilibrium vorticity have a turning point on 
the shell. Some of the rotation laws examined in Paper I 
meet this condition when the degree of differential rotation 
exceeds a certain threshold value. All modes crossing the 
corotation band boundary above this threshold become dy- 
namically unstable. In addition, we observe for one rotation 
law the development of dynamical instability via the merger 
of two stable modes outside the corotation band. 

Paper I raised a number of questions. Is it possible to 
derive a necessary and sufficient condition for instability? 
Do the singular solutions of the continuous spectrum have 
physical meaning? And do the zero-step solutions appear 
distinct from the continuous spectrum in time evolutions? 
In this paper we resolve these questions by addressing the 
initial value problem for the shell. By trying to isolate the 
various contributions to the evolution of smooth initial data, 
we hope to clarify the physical role of different parts of the 
spectrum of a differentially rotating system. The results for 
this simple toy model will guide us in future investigations of 
the t hree dimensional problem iWatts. Andersson fc Joned 



2 KEY RESULTS FROM THE MODE 
ANALYSIS 

In Paper I we analysed the normal mode problem for the 
differentially rotating spherical shell. Since we will refer to 
a number of the results from the mode analysis in this pa- 
per, we summarise the key results in this section. Consider 
an incompressible fluid on a differentially rotating spheri- 



cal shell (with radius R). Assuming that the perturbations 
are cx exp(im(p) we have the ^-component of the perturbed 
Euler equation. 



dtSve + imflSvo — 2Q,cos95vip = —dgSP, (1) 

pit 



and the lys-component 



pR sin 6 

where the equilibrium vorticity is 



(2) 



(3) 



The perturbed continuity equation for an incompressible 
fluid is 



(sinOSve) + imSv^ = 0. 



(4) 



This tells us that only toroidal modes are permitted in this 
problem. Defining the stream function (in the standard way) 



U — Svg sin t 



(5) 



we can combine equations Q and l|2J to give the vorticity 
equation 



where 



{idt - mQ)VtU + -T^U = 0, 
sm 



(6) 



(7) 



sm y sm t) 
is the Laplacian on the unit sphere. Changing variables to 
X — cos the vorticity equation becomes 

{idt-mn)\/lU -mQ'U ^0, (8) 

where 

,2 



= (1 - x')di - 2xd^ 



m 



l-x 



2 ' 



(9) 



and we use a prime to denote derivatives with respect to x. 
In terms of x the equilibrium vorticity is 



n = 2xn- (1 x'^)0.'. 



(10) 



Assuming that we are interested in a mode solution behav- 
ing as exp{—iijjt), where u is the frequency, equation (|5J 
becomes 



{u ~ mQ.)VlU - mn'U = 0. 



(11) 



At corotation points (x = Xc), ui = m^ixc) and the equa- 
tions are formally singular. The nature of the solutions in 
the vicinity of the corotation point can be analysed using 
the method of Frobenius. The general solution is found to 
be 



U = f/.cR + a. 



eg -r "-'sing, 



where 



f^reg — ^ ^ f^n (^ Xc) , 



(12) 



(13) 
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and 



Casing — 



anJyX — XcY* \n{x — Xc) 

+ Cn{x — XcY X > Xc 

b„(xc/x — 1)"'^^ ln{xc/x — 1) 



+d„{xc/x - 1)" 



X < Xc. (14) 



We use separate solutions for Using on either side of the 
corotation point to avoid taking the logarithm of a negative 
quantity. The two solutions for [/sing are matched such that 
U is continuous at the corotation point. The first deriva- 
tives are however discontinuous at Xc, with (in general) both 
a logarithmic singularity and a step discontinuity. The size 
of this step discontinuity in U' varies as one traverses the 
frequency band spanned by the continuous spectrum, and 
for certain special frequencies it vanishes. We refer to the 
solutions at these frequencies as zero-step solutions. At the 
zero-step frequencies the Wronskian of the two linearly in- 
dependent solutions to equation (lllll vanishes, despite the 
logarithmic discontinuity in the first derivative. The impor- 
tance of this will become apparent. 

Let us now consider dynamical instabilities. From equa- 
tion JSJ it is possible to derive the necessary (but not suf- 
ficient) condition for instability Cl' = 0. This condition is 
met for some of the rotation laws considered in Paper I 
when the degree of differential rotation exceeds some thresh- 
old value. Dynamical instabilities are then found to develop 
both where modes meet the corotation band boundary, and 
where modes merge outside the corotation band. 



3 EVOLVING INITIAL DATA: GREEN'S 
FUNCTION APPROACH 

In this paper we consider equation Q as an initial-value 
problem. That is, we are interested in the evolution of pre- 
scribed initial data, say. 



U{x,t = 0) = Uo{x). 



(15) 



We first transform the problem into the frequency domain 
via the Laplace transform (we use s —iio in order to allow 
direct comparison with the mode calculations in Paper I). 



U = 



with inverse 



U = 



2tv 



oo + ic 



Ue '^'^ duj. 



(16) 



(17) 



Applying this transform to equation Q we find 

[uj ~mn)VlU -mQ!U = iVlUQ. (18) 

Note that V^f/o is always non-zero for this problem. The 
only non-zero smooth initial data for which V^I7o — would 
violate the boundary conditions. We should also make a brief 
comment on the Laplace transform. For growing oscillatory 
U, with time dependence exp[(a -I- i/3)t], equation 11811 is 



only correct if the imaginary part of tj, tj/, is greater than 
a. Otherwise the transform would fail to converge as t oo. 
We must therefore ensure that the constant c, which appears 
in the limits of integration of the inverse Laplace transform 
in equation 11711 . lies above all of the poles of the integrand. 
This means that we must have c > a. Thus on the line in- 
tegral, we always have uji > a, and we never breach the as- 
sumption that we made when taking the Laplace transform. 
We then assume analytic continuation of the Laplace trans- 
form into the entire complex uj plane, and use the residue 
theorem to evaluate the line integral. This does not change 
the situation: the line integral itself always has lji > a and 
the Laplace transform and its inverse are well-defined. 

We will now use the Green's function method to derive 
a formal solution to this equation (assuming for the moment 
that the problem is regular). Given a solution G{uj,x,y) to 



G = 5{x-y), 



we have 



U(uj,x) 



1 G{uj,x,y) ^2 
- mQ.[y) ^ 



^lUo{y)dy. 



(19) 



(20) 



Note that because of the symmetry of the problem we need 
only consider the upper hemisphere of the shell. Finally, in- 
verting the transform we obtain the solution to the initial 
value problem 

G{uj,x,y) 



Uix,t) = - 



duo e 



^ - m^l{y) 



y;Uo{y)dy. 

(21) 

The analytic properties of the integrand determine the spec- 
trum of the problem. 

In the standard way, we construct the required Green's 
function G{uj,x, y) from solutions to the homogeneous prob- 
lem. Thus we want to solve 



{to - mQ.)\llU - mtl'tr = 0. 



(22) 



Consider two solutions, Ul and Ur, that satisfy the bound- 
ary conditions at the left and right endpoints of the interval, 
respectively. Using the differential equation one can easily 
show that the Wronskian of any two solutions can be writ- 
ten 



W = UlU'r- UrU'l 



W{uj) 

1 - x2 



(23) 



Pick an arbitrary point y € [0, 1]. Then an overall solution, 
that is continuous ak, x = y, can be written 



U^G 



[h{x - y)UL{y)lJR{x) + H{y - x)U R{y)U l{x)\ 



(24) 

Working out the derivatives of this function, and using equa- 
tion 1221 1. we find that W is a solution to 



■U = 5{x~y)GW{uj). 



(25) 



Comparing this to the equation for the Green's function in 
the previous section we see that we can identify 

Gic.,x,y)^^^^^^ = ^\H{x^y)Ur,{y)Un{x) 
GW{lu) W{uj) L 

+H(y-x)URiy)UL(x)], (26) 
(clearly independent of the normalisation of Ul and Ur). 
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4 SOLVING THE INITIAL VALUE PROBLEM 

To solve the initial value problem analytically we would need 
to evaluate the integrals: 



J{x,y)e 
CW 



[uj - mn{y)] 



-duj 



dy, 



(27) 



(28) 



where U is given by equation (|23J. Unfortunately, we can- 
not in general solve even the homogeneous problem (for Ul 
and Ur) analytically. It is however possible to generate spe- 
cific toy problems for which the initial value problem may 
be solved analytically. We are particularly interested in the 
behaviour of the continuous spectrum and the zero-step so- 
lutions. We therefore develop two toy problems. The first 
is a particular rotation law that admits only a continuous 
spectrum as a solution. An analytic solution of the initial 
value problem is possible for simple initial data, and the 
solution sheds light on the physical nature of the perturba- 
tion. In the second problem we aim to investigate the nature 
of the zero-step solutions. Ideally one would like to find a 
simple rotation law, with a zero-step solution, for which one 
could solve the initial value problem analytically. We have 
not been able to identify such a law. Instead we have de- 
veloped a toy mathematical problem with solutions that in- 
clude both a continuous spectrum and a zero-step solution. 
The behaviour of the eigenfunctions of this system near the 
critical point is mathematically identical to the behaviour 
that we find for the shell problem. The initial value prob- 
lem for our toy system can however be solved analytically, 
guiding us towards a quasi-analytic solution for the zero-step 
solutions and continuous spectrum of the shell problem. The 
details of the two toy problems are contained in Appendices 
^and^ In this section we summarise only the key results 
from each problem. 

Appendix ^presents a full analytic solution of the ini- 
tial value problem on the shell, for a simple rotation law 
that possesses only a continuous spectrum (with no zero- 
step solutions). For simple initial data the collective pertur- 
bation associated with the continuous spectrum is found to 
be non-singular, and hence physical. The perturbation dies 
away with time, with inverse power law decay rather than 
exponential decay. The frequency of the decaying perturba- 
tion depends in part upon position. Similar behaviour has 
been found for t he continuous sp ectrum of a differentially 
rotating cylinder (lBalbinskilll984f) . 

Appendix^introduces a toy problem that has singular 
continuous spectrum eigenfunctions with behaviour similar 
to those of our problem. By careful selection of boundary 
conditions, one can ensure the presence of a single zero-step 
solution within the continuous spectrum. In this case we are 
able to make some progress towards an analytic solution of 
the initial value problem. We draw several important conclu- 
sions that are likely to transfer to the differentially rotating 
shell. Firstly, the physical perturbation contains a constant 
amplitude oscillatory component at the zero-step frequency. 
Secondly, just as in Appendix^ the rest of the continuous 
spectrum dies away with time for simple initial data. Again 
we find inverse power law decay rather than exponential de- 



cay. The frequency of the perturbation associated with the 
continuous spectrum is, as before, found to depend in part 
upon position. 

Let us now proceed as far as we can with the analytic so- 
lution to the initial value problem for the differentially rotat- 
ing shell. We start by considering the integral over position, 
equation 1281 . We will write uj — mQ{y) = Q.{uj,y){y — yc), 
where the corotation point t/c is a function of lu. Equation 
12811 becomes 



J{x,y) 



Ul{x) 



+Ur{x) 



1 UR{y)\/lUoiy) 



^{^,y){y-yc) 
ULiy)VlUo{y) 



dy 



dy. 



(29) 



/o ^{^,y){y ~ yc) 

If u! is in the range of niQ, then < j/c < 1 and the integrand 
will be singular a.t y = yc. In Paper I we examined the 
variation of U with y and found that in the vicinity of the 
corotation point {y « yc), for real ui, it behaves as 

oo oo 

U = Y1 Cn{y - yc)" + a„{y - yc)"^^ ln(i/ - yc). (30) 

Whether yc is greater than or less than x will depend on 
Lu. If J/c > X, the first integral in equation 1291 has a pole 
and a branch point at y = yc, and we must treat the in- 
tegral as a principal value integral. If on the other hand 
yc < X, then it is the second integral that must be treated 
as a principal value integral. The same procedure is followed 
in Appendix|B| These integrals are not easily evaluated ana- 
lytically; one must in general resort to numerical methods to 
solve even the homogeneous problem for Ul and Ur. How- 
ever, by examining the power series expansions around the 
singular point, which are identical to the power series expan- 
sions around the singular point for the problem in Appendix 
[E\ it is easy to see that the principal value integrals must 
result in smooth, non-singular functions. Thus even with- 
out solving the integrals analytically, we can state that the 
first integral will give rise to a function Ii{x,uj) and the sec- 
ond integral to a function I2{x,lj), where both 7i and I2 
are smooth, non-singular functions. The frequency integral, 
equation 1271 . becomes 



U{x,t) = - 



2n 



^c W{u) 



\—\UL{uj,x)h{x,U)) 



+ UR{L0,x)l2{x,L0)\e '"^ diU 



(31) 



We cannot solve this integral analytically. Instead, we will 
use Cauchy's theorem to extract its properties. We complete 
the contour using a semicircle of infinite radius in the lower 
half plane that must be deformed around any branch cuts. 
The contribution from the semicircle integral vanishes be- 
cause the integrand tends to zero on this part of the contour. 
The only non-zero contribution will come from the poles and 
branch cuts of the integrand. 

We know from Paper I that there are zeroes of W at 
frequencies at which there are discrete normal modes, and 
at frequencies within the continuous spectrum where there 
are zero-step solutions. W has a zero of multiplicity one in 
all cases except where modes or zero-step solutions merge: 
at these frequencies it has a zero of multiplicity two. 
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We must also consider the behaviour of U . This function 
is well-behaved and non-singular except for uj within the 
range of mf2. Let us first consider the behaviour near the 
frequency uj = mQ,[x), the corotation frequency. We need to 
know how U varies with uo for fixed x. This would in principle 
require a numerical calculation, but for the fact that we do 
know how U varies with x for fixed to in the vicinity of 
the point x = Xc- Because Xc is a function oi uj, we can 
rewrite our power series expansions as an expansion around 
the singular frequency rather than the singular position (see 
Appendix|n| where the equivalence is clearer) . We find that 
U is non-singular but has a branch point at a; = niQ,{x). 
The two end points of the continuous spectrum, uj = mQ.Ynm 
and UJ = mJlmax are also branch points. The contour of 
integration should look similar to that shown in Figure IbTI 

We can now make some definitive statements about the 
time-dependence of the physical perturbation U{x,t). 

(i) The poles of the integrand at discrete frequencies 
where the Wronskian is zero will contribute terms with sim- 
ple oscillatory time dependence exp{iujrit). This is true both 
for discrete normal modes and the zero-step solutions. Thus 
by virtue of having zero Wronskian, the zero-step frequen- 
cies should appear virtually indistinguishable from the stan- 
dard discrete normal modes, despite being in the continuous 
spectrum. 

(ii) At frequencies where modes merge and the Wronskian 
has a double pole we will have additional terms of the form 
texp{iujnt). Such linear growth signals the onset of instabil- 
ity, as can be seen if one considers an unstable mode ex- 
hibiting exponential growth exp[(itt;„ +ujm)t]. For small 
this can be expressed as a Taylor series: 



We assume that A is real (note that, strictly speaking, A 
will depend on uj and y). Then we find from equation 1251 
that 

A^CW+ r ""^' udx. (34) 
Jq uj — mil 

Hence 

€ = CW = A- / ^ Udx. (35) 

Jo UJ — mil 

In the limit when uj is real but outside the range m^l{xc) 
for Xc in [0,1], e is real, e has zeroes at frequencies where there 
are normal mode solutions outside the corotation band. If 
on the other hand uj is real, but in the range uj = mQ,{xc) 
for Xc G [0, 1], then the integrand is singular and the integral 
must be interpreted as the principal value 

e = CW = A-P f '^^ ^ U(uj,x,y)dx 
Jq uj~m\l 

±«7r^^7|^W(a;,Xc,2/). (36) 

The sign depends on whether we circle the singularity in the 
upper or the lower half-plane. This quantity e, as defined in 
equation (I36II . can be used to derive a sufficient instability 
criterion. Our next step is to derive a useful expression for 
U. 

In Section @ we constructed a Green's function for the 
initial value problem, equation 12611 . that depended on U, 
which depends in turn on the solutions, U , to the homoge- 
neous problem. The equation for the homogeneous problem, 
equation 1221 may be rewritten as 



'[1+Umt+ ...]. 



(32) 



This behaviour is to be expect ed since mod e mergers usually 
lead to dynamical instability JSchutdllQSd) . 

(iii) The behaviour of the continuous spectrum is more 
difficult to ascertain but, since the end points of the corota- 
tion band on the real line are branch points, we would expect 
decaying oscillatory behaviour similar to that found in Ap- 
pendix^ The frequency of oscillation will inevitably depend 
in part upon position. We cannot rule out the possibility of 
constant amplitude oscillation, or even growth of the con- 
tinuous spectrum, in response to some initial data. What is 
clear, however, is that the perturbation associated with the 
continuous spectrum is non-singular and hence physical. 



UJ — mil 



(37) 



We can solve this equation using the Green's function 
Go{x,x') that satisfies 



V^Go = 5{x — x). 



(38) 



The explicit form for this Green's function is given in Ap- 
pendix^ Finally, we have 



U{x) 



^ mO.' 
UJ — mQ 



U Go{x, x')dx' 



(39) 



Knowing that we can find Go , we can rewrite equation 1251 
for U as an integral equation. We obtain 



5 A SUFFICIENT INSTABILITY CONDITION 

Having discussed the general properties of the solution, we 
will now use an analytic approach to the problem to derive 
a sufficient condition for instability. We follow the method 
outlined bv lBalmforth fc Morrison! ^1999^ for inviscid plane 
parallel shear fiow. We start by reference to equation 12511 . 
the equation obeyed by W. Let us assume that we can nor- 
malise the solution in such a way that 

/ VlUdx = A. (33) 
Jo 



U{uj,x,y) = CWGo{x,y) 

+ r ^^^^U{Lo,x',y)Go{x,x')dx', (40) 
Jq uj — mil 

or, if we use the normalisation, equation 1331 , 

^ If the eigenfunction U were regular at the corotation point then 
the integrand would not be singular. This applies to both types of 
regular eigenfunction described in Paper I, and can be seen from 
the power series expansions of the regular eigenfunctions around 
the singular point. 
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l^{^,x,y) = AGo{x,y)+ [ r, ^{^,x',y) 

Jq lo — mil 



[Go{x,x ) — Go{x,y)]dx . 



(41) 



Recall that we have so far assumed that uj is complex in 
order to avoid singularities on the real a;-axis. If we relax 
this assumption, we need to allow for the possibility that 
UJ — m^{xc) = 0. Then we find that I41II is a singular prob- 
lem. However, we can regularise this equation by a judicious 
choice of y. We therefore take y = Xc in 14111 to get 



U{lJ,X,Xc) 



AGo{x,Xc) + 



-U {u, X , Xc 



[Go{x,x')-Go{x,Xc)]dx'. (42) 

As discussed by iBalmforth fc MorrisonI 1^9^ , a solution 
to the Fredholm equation (|42^ is a singular eigenfunction 
associated with a continuous spectrum frequency. Note that 
we can let U{uj, x, Xc) — > U{uj,x) since Xc is determined by 
UJ anyway. 

Equations 1361 and 14211 hold for all points where the 
normalisation given by equation 1331 is valid. There are how- 
ever certain solutions within the continuous spectrum for 
which this normalisation is not appropriate. If we have solu- 
tions for which W = Q then equation 12511 for U is identical 
to the normal mode equation, Ijll^ . that we solved in Pa- 
per I. Thus for zero Wronskian solutions we can identify U 
with the normal mode solutions found in Paper I. Most zero 
Wronskian solutions correspond to modes with finite and 
continuous eigenfunctions outside the corotation band. The 
zero-step solutions are an exception, having zero Wronskian 
and a logarithmic discontinuity in the first derivative. So for 
the frequencies where zero-step solutions exist we must treat 
the integral in equation 1331 as a principal value, with 



A^/ VlUdx = V ( VlUdx±iTi[{x- Xc)Vlu\ 
Jo Jo 



For these special frequencies, 



(43) 



r 



-u 



[{x^x,)Vlu]^ + 



dx 



fi'fXe) 



U{uj,xc,y) 



(44) 



That the real part of e should be zero if W is a solution to the 
normal mode problem is clear from equation Ull^ . That the 
imaginary part must also be zero can be seen by taking the 
limit aax ^ Xc of {x — Xc) times equation IIH . and applying 
L'Hopital's rule. So as we would expect, e (and hence W) 
is zero at frequencies where there are zero-step solutions to 
the normal mode problem. For all other frequencies within 
the continuous spectrum e will in general be complex and 
non-zero. 

At frequencies where there are zero-step solutions, the 
regularising term in the integral in equation 1421 cancels 
exactly the term involving A, and equation 1421 becomes 



Uiu^x) 



UJ — mQ 



U{uj,x')Go{x, x')dx' . 



(45) 



The equation becomes homogeneous, with a kernel that is 
truly singular. The presence of homogeneous solutions at 
certain frequencies complicates any attempts to solve equa- 
tion 1421 directly. Fortunately, we do not need to solve this 
equation to derive a sufficient instability criterion. 

Let us recap. Given a solution to the integral equation 
1421 . we can find the quantity e defined in equation 1361 . 
These two equations hold for all frequencies in the continu- 
ous spectrum apart from those for which there are zero-step 
continuous spectrum solutions to the normal mode problem. 
At these frequencies, e will be zero and there will be a ho- 
mogeneous solution to equation 1421 . as discussed above. 

Let us rescale equations 1421 and 1361 according to 
U/A U and e/A ^ e to get 



U{uj,x) = Go(x,Xc)+ [ ^ U(uj,x') 
Jo UJ — mil 



and 



r 



UJ — mil 
[Go{x,x') — Go{x,Xc)]dx' , 



mO! , , , , , Q.' (xc) ^ ,, 

tU(uj, x)dx ± iT^TTT-f — rM(cj, Xc 



(46) 



(47) 



n'{xc) 

In order to have an instability we must have a zero of the 
Wronskian W in the upper half of the oj-plane. The Nyquist 
method tells us that the number of zeroes A/" of W in a 
region enclosed by a contour C in the uj plane is given by 
the integral 



(48) 



where A is used to denote the change in a quantity. The 
integral in equation 1)48^ tells us the number of times that 
the function W winds around the origin of the complex W 
plane as uj varies around C. We decompose the integral into 
two parts: 



1 f W\ 1 f W , 1 f W , 

/ —duj = / —duj H / —duj. (49) 

2ni Jc W 27ri J^ W 2ni Jc, W 

The first integral is taken over the real line TZ, the second 
over a semicircle C' in the upper half plane. In Appendix Icl 
we show that the integral over the semicircle vanishes. Now 
e{oj) = G{uj)W{uj), so 

W' f e' f C' 

-^^duj — / —duj — / —duj. (50) 
m W Jtz ^ Jn C 

C was an arbitrary constant in equation 1241 . so we can 
choose C to be real. Then, setting e = |£| exp(iarge) and 
= IW"! exp(iargM^), equation 1481 becomes 



AA= — Aarg(e). 



(51) 



In other words, the number of unstable eigenvalues should 
be determined by the change in the argument of e as lj tra- 
verses the real line. In fact, we can constrain the range of 
integration still further. For frequencies outside the coro- 
tation band, the principal value integral in equation 1471 
becomes a regular integral with no imaginary part (since we 
can normalise U to ensure that it is real). Hence, e is real 
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for all frequencies outside the range mf^min < u) < mf^max 
and cannot wind. The presence of a corotation band and 
the associated continuous spectrum is clearly essential for 
instability on the shell. Thus 



AA= — Aarg(e) = 



— duo. 



(52) 



The existence of e = solutions at certain frequencies 
within the corotation band, and for certain degrees of dif- 
ferential rotation on the edges of the corotation band, will 
not affect the winding number. This is because these are iso- 
lated points where e jumps directly to zero; and there is no 
winding if e passes through the origin rather than around it. 
The frequencies where e = can be excluded from the ca l- 
culation of winding number jBalmforth fc Morrisonlll99El) . 
and will not be discussed further. 

At the endpoints of the corotation band, lu = mf2min 
and UJ = mfimax, one can demonstrate that e = 1. This 
means that there can only be winding if we cross the negative 
real e-axis while traversing the band. 

From equation 14711 . it is clear that the imaginary part 
of e is zero if either U{u,Xc) = or Cl'{xc) = 0. The first 
condition, U(ijj, Xc) = 0, cannot lead to instability, because if 
this is the case we have a regular mode within the corotation 
band, and the real part of e would also vanish. Thus for 
winding we require 



q!{x^) = 0, 



(53) 



which is consistent with the necessary criterion for instabil- 
ity derived in Section |S] We denote the frequency at which 
this inflexion point condition is met as uji. 

For there to be at least one unstable eigenvalue, the real 
part of e must be negative at the frequency Ui . We therefore 
have a sufficient condition for instability 



V 



mO.' 



U{ijJi,x)dx > 1. 



(54) 



The above condition tells us that there is at least one un- 
stable eigenvalue. But what if we want to know the absolute 
number of unstable eigenvalues? In Paper I we found the 
total number of unstable eigenvalues by computing normal 
modes. The method outlined in this section offers in prin- 
ciple an alternative. If we could solve equation (I46II numer- 
ically, using a Fredholm equation solving routine, we could 
use the resulting U to compute e as the frequency varies 
across the continuous spectrum range, and hence determine 
the winding n umber directly. A similar calcu lation has been 
carried out bv lBalmforth fc MorrisonI il999l) for plane par- 
allel flow. Their problem is however simpler as there are 
no homogeneous solutions to their equivalent of equation 
14611 . In our case one would have to begin by solving the 
homogeneous (eigenvalue) Fredholm problem, as stated in 
equation 1451 1. and then exclude these frequencies from the 
main calculation. We have assessed the feasibility of solving 
the homogeneous equation using the n umerical methods for 
singul ar integral equations discussed in lDelves fc Mohamedl 
il985h . but conclude that application of these methods is 
far from trivial in this case. Computation of normal modes 



is a far more convenient method of determining the total 
number of unstable eigenvalues for this problem. 



6 EVOLVING INITIAL DATA: NUMERICAL 
RESULTS 

The prime way to check the mode calculations is the numer- 
ical evolution of generic smooth initial data. This will also 
enable us to check the quasi-analytic solution to the initial 
value problem formulated in Section 3] 

We start by reformulating equation Q to ensure that 
our variables are real valued. We introduce 

U{t,e,(p) ^ C{t,9) cos mip- S{t,e) sin mip (55) 

After decoupling the sine/cosine parts, we arrive at two 
equations that can be written succinctly as 



dtA - mflB + mdeflS = 
dtB + mflA - mdedC = 

where A and B are defined by 

Vie = ^ 
sm 6 

Vis = 

sm 6 

In order to get a first-order system we introduce 

Z = deC 
X = dgS 

so equations 15811 and I59I I become 
sin — sin^ 6*98^ — sin 6* cos -f m^C — 



sin 6B — sin OdeX — sin 6 cos 9X + m S 



(56) 
(57) 



(58) 
(59) 

(60) 
(61) 

(62) 
(63) 



We now have six first order equations for six unknowns. We 
can see from equations H62^ - I|63|l that for regularity at the 
poles we require C — S = 0. 

Our numerical solution of the six equations is based on 
a second order accurate Crank-Nicholson implicit finite dif- 
ferencing scheme. The finite differenced equations are solved 
using a relaxation scheme with 601 gridpoints and a timestep 
At — 0.05AS. Convergence testing has verified that the code 
is second order accurate, and stability analysis indicates that 
the code is stable over many dynamical timescales. To test 
the code we ran time evolutions for the simple rotation law 
described in Appendix^ The analytical solution of the ini- 
tial value problem was used to verify the robustness of the 
numerical time evolution code. We found the results to agree 
extremely well with those shown in Figure IXTl 

We then performed time evolutions for a more compli- 
cated rotation law, the solar rotation law that we studied in 
Paper I (Law 3 of that paper): 

n = 27r (454.8 - 60.4/3 cos^ 6* -71.4/3 cos* nHz) (64) 

The parameter /3 measures the degree of differential rota- 
tion, uniform rotation corresponding to /3 = 0. 

For initial data we used the standard spherical harmon- 
ics Yi"^ that are the eigenfunctions of U in uniform rotation. 
A sample of numerical results for this rotation law with 
m = 2 are shown in Figures [1121 For frequencies outside the 
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f(nHz) f(nH2) f(nHz) f(nHz) 



Figure 1. Power spectrum of numerical time evolution of data for 
the solar differential rotation law with m = 2 and f} = 0.4. The 
four panels represent initial data proportional to the spherical 
harmonic , with Z = 2 — 5. The sampling position is 9 = 7r/30. 
The dashed line marks the lower edge of the corotation band, 
at 804.6. We know from Paper I that at /3 = 0.4 the I = 2, 
3, and 4 modes are outside the corotation band, at frequencies 
592.1, 744.3 and 796.0 respectively. Peaks at these frequencies are 
clearly visible. The Z = 5 mode does not cross into the corotation 
band but the mode results predict a zero-step solution just inside 
corotation, at frequency 806.35. A peak at this point is visible in 
the power spectrum. 



Figure 2. Power spectrum of numerical time evolution of data for 
the solar differential rotation law with m = 2 and /3 = 0.5. The 
four panels represent initial data proportional to the spherical 
harmonic , with Z = 2 — 5. The sampling position is 9 = 7r/30. 
The dashed line marks the lower edge of the corotation band, at 
777.8. We know from Paper I that at /3 = 0.5 the / = 2 and 3 
modes are outside the corotation band, at frequencies 588.66 and 
741.08 respectively. Peaks at these frequencies are clearly visible. 
The I = 4 mode is now inside the corotation band, at frequency 
788.34. This is visible, in addition to a zero-step solution excited 
by / = 5 initial data at 802.50. The peak at around 825 for I = 
4 and 5 initial data is probably associated with the continuous 
spectrum and is an artefact of sampling position. 



corotation band, there are clear peaks at the frequencies of 
the regular normal modes. For dynamically unstable modes 
we are able to verify the growth times predicted by the nor- 
mal mode analysis. What is most interesting, however, is 
the behaviour of the continuous spectrum. The power spec- 
trum shows distinct peaks within the continuous spectrum 
at frequencies where we found zero-step solutions in Paper 
I, just as predicted in Section |1] Although the size of the 
peak in the power spectrum varies depending on sampling 
position (the value of x at which we monitor the time evo- 
lution), these frequencies clearly stand out from the rest of 
the continuous spectrum. 

We were also interested to see how the rest of the con- 
tiimous spectrum behaved in response to an initial pertur- 
bation. The continuous spectrum was observed to die away 
with time for smooth initial data, and to be excited to a 
much lower extent than the modes and zero-step solutions. 

The numerical evolutions support the conclusions from 
our Green's function analysis. In particular, it is clear that 
the zero-step solutions stand out from the rest of the con- 
tinuous spectrum. Hence, it is reasonable to expect them to 
be of particular relevance for the dynamic al evolution. Thi s 
would be in accordance with the results of lBalbinskil il985l) , 
which show that the merger of zero-step solutions may lead 
to dynamical instabilities. 



7 DISCUSSION 

The main question arising from the normal mode analysis of 
Paper I was the nature of the physical perturbations associ- 
ated with the singular eigenfunctions of both the zero-step 
corotating solutions, and the continuous spectrum within 



which they exist. By studying the initial value problem both 
analytically and numerically we have been able to resolve 
this issue. 

Firstly, we have confirmed that the physical pertur- 
bations associated with these phenomena are not singular. 
Their time dependence in response to an initial perturba- 
tion, however, is complicated. For appropriate initial data, 
the zero-step solutions behave in much the same way as sta- 
ble normal modes outside corotation. That is, they produce 
a clear peak in the power spectrum at a fixed frequency, 
standing out from the rest of the continuous spectrum. The 
fact that such oscillatory solutions can exist within the con- 
tinuous spectrum is of great relevance for future studies of 
the oscillations of differentially rotating systems. 

The behaviour of the rest of the continuous spectrum 
is difficult to predict analytically, and can only be eluci- 
dated by numerical time evolutions. The analytic calcula- 
tions suggest that the physical perturbation associated with 
the continuous spectrum is oscillatory, with a complicated 
frequency dependence that varies in part according to the 
position on the shell. The amplitude of the oscillations could 
grow, remain at constant amplitude or die away with time 
(as a an inverse power of time), depending on the initial data 
and the rotation law. For the rotation laws and initial data 
that we have studied we observe only decay, but we cannot 
rule out either constant amplitudes or growth. 

The results of numerical time evolutions have verified 
to high accuracy the predictions of mode frequencies made 
by the normal mode analysis. This is true both for modes 
outside the corotation band (including dynamically unsta- 
ble modes), and the zero-step solutions within it. This is 
encouraging, as when we move on to study more realistic 
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stellar models, we are likely to rely more on time evolutions 
than on increasingly complicated mode analysis. 



the response of the continuous spectrum to an initial per- 
turbation. We select the simple rotation law 
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1 



(Al) 



for which Q! = Q. Solving the normal mode problem for this 
rotation law, one finds that the only solutions form a con- 
tinuous spectrum with derivatives that are discontinuous at 
the corotation point. In contrast to the general case, how- 
ever, the derivatives possess only a step discontinuity at the 
corotation point (no logarithmic singularity). There are no 
zero-step solutions. 

Returning to the initial value problem, equation 1181 
becomes 
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APPENDIX A: A ROTATION LAW THAT 
EXHIBITS ONLY A CONTINUOUS SPECTRUM 

Solving the initial value problem analytically for an arbi- 
trary rotation law is extremely difficult. For certain rota- 
tion laws, however, an analytical solution of the initial value 
problem is both feasible and instructive. In this Appendix 
we analyse a rotation law that exhibits a continuous spec- 
trum but no other type of oscillation. Our aim is to isolate 



(A2) 



We solve this equation using a Green's function Go{^,y) 
that satisfies 



(A3) 



Because of the symmetry of the governing equation there 
are two possible boundary conditions at a; = 0, either 
dGo/dx{Q,y) — (even solutions), or Go(0, j/) — (odd 
solutions). For even solutions, we find that 



Go(x,y) 



and for odd solutions. 



2m 






1 - 

2m 




(f+i) 



Goix,y) 



1 

2m 



1 
2m 



y+1 \ 2 
a-1 



x-l ] 2 _ / x + lV 
x+1 J 



x-1 \ 2 
x + 1 J 

y-l\ 



X > y 

X <y, 
(A4) 



X > y 

X <y. 
(A5) 



We can now recover the physical perturbation using the in- 
version integral 



U{x,t) = / n{x,y)V'yUo(y)dy, 
Jo 



where 



H{x,y) = 



_1_ 

2^ 



Go(uj,x,y)e " 
LO — mQ{y) 



-du. 



(A6) 



(A7) 



We illustrate the typical behaviour by giving results for the 
particular case Uo{y) — Pi^(y), with I — m (even). In this 
case, VlUoiy) oc (1 - j/2)'"/2 ^nd we find 



U{x,t) = 



>m+l 



1 - X\ 2 



47rm I V 1 + 3; 
(l-y)™|e-*'""(''>*dy-F 
1-x 



[{i + yY' 



l+x 



+ 



+ 



l+x 



1 - x^ 

(l-2/)'"e-"""(^''dyl. (A8) 
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Evaluating the integral analytically for Z = m = 1 we find 



U{x,t) 



where 



2b 



g-»n<:t| {a + b) 



(4 + 2it)e-'^ + if - 4if)Ei 



+ 



i - 3a + fe + 2(& - a)a; - i{a + 6)(f + x)t 



+ {a + h)x^\e TFF + 
\[A^at-[a + h)t-] Ei (i^ 



(A9) 




Figure Al. Time evolution of \U\ for I = m = 1 and I = 2,m = 1 
initial data. At late times, \U\ falls off as f/t for the I = m = I 
data and as f/t^ for the i = 2,m = f data. 



i(x) = 



1 + a; 



1+x 



(AlO) 



and the asymptotic expansion of the exponential integral 
Ei(a;) is 



Ei(2;) 



(All) 



Expanding the exponential integrals we find that all terms 
whose amplitudes grow with t cancel. The resulting time 
dependence is illustrated in Figure IXTl The most persistent 
term has an amplitude that decays as 1/t. In addition, the 
three different frequency components Qc, O.c/2 and Q,c/{'i + 
x) interact to produce beating. The beating is particularly 
pronounced if x is chosen such that f2c/(l + x) is very close 
to either of the other two frequencies. 

If instead we consider initial data Uo{y) ~ Pi^{y) with 
l = m + l (odd) then V^C/o(?y) oc y{l - y2)™/2 fi^d 
that 



U{x,t) = 



AzTTm I \ I + x 
(l + J/)'"|e-"""(«'*dy + 



1 - X 
1 + x 



1±| j ' £ y{l - j/)'"e-"""(«)*dy}(A12) 

Evaluating the integral analytically, we obtain similar results 
to those for / — m initial data, except that now the most 
persistent terms decay as l/t^. This is illustrated in Figure 

mi 



in addition to exhibiting both of these features, is simple 
enough to permit analytical solution of the initial value 
problem. In this Appendix we instead analyse a toy problem 
that has solutions with the required behaviour, to inform our 
investigation of the shell problem. 



Bl Selection of a toy problem 

Consider the equation 



{x~idt)U"{x) + U{x) = 0, 



(Bl) 



where primes indicate derivatives with respect to x. Assume 
a normal mode solution, with U oc exp{—iujt). Equation 
lIBll l becomes 



{x - uj)U" + U = 0. 



(B2) 



The critical point x = lu, which gives rise to the continuous 
spectrum, plays the same role as the corotation point on the 
differentially rotating shell. The general solution to equation 
llE2l is 



u ^ v^r^[AiJi{2v^r^) + A2Yi{2v^r^)], (bs) 

where and are Bessel functions of the first and sec- 
ond kind respectively, u being the index of the function. Ai 
and A2 are constants. To show that the solutions have the 
required behaviour, it is convenient to consider separate so- 
lutions on either side of the critical point (compare to Paper 
I): 



APPENDIX B: A TOY PROBLEM WITH A 
CONTINUOUS SPECTRUM AND A 
ZERO-STEP SOLUTION 

In Appendix^ we solved the initial value problem for the 
differentially rotating shell for a rotation law that exhibits 
only a continuous spectrum. The continuous spectrum eigen- 
functions obtained for this rotation law differ from those 
found for more general rotation laws, however, in two ways. 
Firstly, the eigenfunctions do not possess a logarithmic sin- 
gularity in the first derivatives, only a finite step. Secondly, 
there are no zero-step solutions. 

We have not been able to identify a rotation law that. 



uix)^l f'<;^; (B4) 

^ ^ \ U>{x) x>iu, ^ ' 

where 

U<{x) = V^^[ciJi(2v/;r^) +C2yi(2V^^)](B5) 
U>{x) = Var^[c3Ji(2V^r^) +C4yi(2V^r^)](B6) 
where ci - C4 are constants. Near the critical point: 

[/< = -— + 0{x~uj) (B7) 

IT 

U> = -— + 0{x-uj). (B8) 
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For continuity of the function at the critical point we re- 
quire C2 = C4. Let us now consider the behaviour of the first 
derivatives near the critical point: 



[/< = ci + —[ln\x-iu\+in + 2'y]+0{x-iL}) 



(B9) 
(BIO) 



U'y = C3 + —[ln\x- 1^1+2-/]+ 0{x-uj), 
n 

where 7 is Euler's constant. Thus at the critical point there 
will be both a logarithmic singularity and a finite step in 
the first derivative. The degree of freedom associated with 
the finite step in the first derivative allows us to construct 
a continuous spectrum of solutions for u) in the range of 
X. Mathematically the behaviour is identical to that of the 
continuous spectrum on the differentially rotating shell. The 
Wronskian of the two solutions, W, is 



We will solve this problem using a Green's function that 
obeys the equation 



G" + 



G 



S{x - y). 



(B19) 



B3 Finding the Green's function 

We construct the Green's function as follows. Let Ul and 
Ur be solutions to the homogeneous equation that obey the 
boundary conditions on the left and right of the domain 
respectively. The Wronskian of these two solutions, W , is a 
function of frequency only. Consider the function 



W = — [C4C1 — C2C3 +iC2Ci\. 
TT 



(Bll) 



Since C2 = C4 we see that W is zero if ci — C3 + ic2 — 0. 
This corresponds to the finite step in the first derivative 
vanishing, indicating the presence of a zero-step solution. 
As we are interested in the behaviour of zero-step solutions 
we will select a domain and boundary conditions such that 
there is at least one zero-step solution. We select the domain 
X = [0, 1] and fix one of the boundary conditions, ?7< (0) — 0. 
In this case there is a zero-step solution for uj = 1/2 if we 
insist on a second boundary condition (7>(1) ~ P, where 



C3 lim 



-[/2/3-/l/4 + i/l/3] 



/2 - ih 



(B12) 



and 




(B13) 
(B14) 
(B15) 
(B16) 

There are no other zero-step solutions within the continu- 
ous spectrum with < u; < 1. Investigation suggests that 
there are no other mode solutions for this problem with these 
boundary conditions with either complex uj or real u) outside 
the range < a; < 1. 



F{x,y) = C[H{y - x)U R{y)tj l{x) + H{x ~ y)U L{y)\J r(x)\. 

(B20) 

F is continuous at the point x = y, irrespective of the nor- 
malisation of Ol and Ur. It can be shown to obey the equa- 
tion 



= CW{Lo)S{x-y). 



X — LO 

We can see that a valid solution to equation l|B19^ is 

G{u,x,y) = ^^^[H{y-x)UR{y)UL{x) 
+H(x-y)UL{y)UR(x) 



(B21) 



(B22) 



We now need to determine Ul and Ur. From our examina- 
tion of the normal mode problem, we know that normalised 
solutions to the homogeneous problem are 

Ul{x) — ^/x~~uj [Ji{2Vx^) + C2Yi(2V;F^)] (B23) 



Ur{x) = ^/^^ [Ji(2v/^~;:J) -I- C4yi(2x/^^)] , (B24) 

where (comparing to the normal mode problem) we have de- 
fined C2 = C2/C1, C4 = C4/C3, and /3 = /3/c3. The Wronskian 
of the two solutions is 



B2 Setting up the initial value problem 

In the previous section we established the boundary condi- 
tions necessary to ensure that the toy problem exhibits both 
a continuous spectrum and a zero-step solution at ci; = 1/2. 
We now define 



'0 



(B17) 



and apply a Laplace transform to equation IBH to obtain 
the inhomogeneous equation 



U 



iUk 



(B18) 



W = Ci - C2 + ic2CA. (B25) 
The boundary condition at a:: = 0, Ul{Q) ~ 0, gives 



C2 = 



Jl(2^/=ZJ) 



Yi{2^) 

The boundary condition at a; = 1, Ur{1) = (3, gives 



Ci = 



Yi(2VT^) 



Ji{2y/T^) 



(B26) 



(B27) 



It can be verified that the Wronskian is indeed zero at uj 
1/2. We have now found our Green's function. 
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B4 The inversion integrals (part I) 

The first stage in recovering U{x,t) is to use tlie Green's 
function to find U{x,lu). Using equations ljB18|l and ljB19fl 
we can sliow, in tlie standard way, tliat 



Uix) 



y-^ 



dy 



[cr/' - 



UG' 



(B28) 



Tlie last term in equation lB28t is non-zero because of our 
choice of boundary conditions. By choosing conditions that 
ensured the presence of a zero-step solution, we have also 
ensured that G and U do not obey the same boundary con- 
ditions. In the following sections we will neglect this second 
term because our primary interest is in the zero-step so- 
lutions. Their contribution is encapsulated in the position 
integral, not in the second term, which depends solely on 
the end-points. Restricting our attention to the position in- 
tegral. 



Ul>'{y)GiLJ,x,y) 
y-uj 



dy 



+ Ur{x) {Ip3 + C2{Lu)Ip4) te29) 



where 
Ipi ~ 

Ip2 = 
Ip3 ~ 

Ipi = 



Ji(2^y — uj) 



Y^^VyZR 
Vy-^ 



dy 


(B30) 


dy 


(B31) 


dy 


(B32) 


dy. 


(B33) 



In order to evaluate these integrals we must select initial 
data that obeys the boundary conditions. We will consider 
the simple case 



[/o(y) = 2/72 + (/9- 1/2)2/ 



(B34) 



Ipi and Ip3 are regular at the point y — lj, so the integrals 
are straightforward: 



Jo{2y/I^) - Joi2y/l -Lo), 



IpS = -Jo{2Vx-Lo) + Jo(2v^). 



(B35) 



(B36) 



The integrals Ip2 and Ipi require a little more care because 
the integrands are singular at the point y = uj (real uj only). 
\i X < UJ < 1, we must treat Ip2 as a principal value integral, 
and Ipi is non-singular. If on the other hand Q < uj < x, 
we must treat Ipi as a principal value integral, and Ip2 is 
non-singular. Thus we get 



lp2 



Ipi 



yo(2V^r^) - Yo{2^^T^) - H{uj - x)i, (B37) 



-Yo{2V^^^)+Yo{2./^)- H{x-uj)i. (B38) 




Figure Bl. The contour for the frequency inversion integral. 
Filled dots indicate branch points, open dots isolated poles. The 
thick solid lines are branch cuts. 



B5 The inversion integrals (part II) 

The final inversion integral to recover the time dependent 
solution U(x,t) is 



U{x,t) 



1 

2^ 



U{uj,x)e '"'dLj. 



(B39) 



For our given choice of initial data, equation IIB39II becomes: 



U(x,t) 



1 

2^ 



oo+ic 



■ I Ji {2Vx^) 



oo+^c W{UJ) 

Jo(2V^) - Jo(2yr^)] + Ji(2y X — UJ) 



[(C4 - C2)yi,(2V^^) - 04^0(2^1^) - 

C2Yo{2^)] + yi(2V^^)[(c2 - Ci) 

jo{2V^r^) + ciJo{2y/~^) - c2Jo(2yr 

—ic2C4Y-i_{2\/lE~~uj) — iJ\(2\Jx — uj) 

\II{UJ - x)Ci + II(X - UJ)C2\ |. 



(B40) 



The integrand has a simple pole at = 1/2 due to the 
zero of the Wronskian. In addition, there are branch points 
at ci; = Ojj:,!. The contour of integration for x < 1/2 is 
shown in Figure IHTl fif x > 1/2, we would instead join the 
points u) = X and u = 1 by a branch cut, and have an 
isolated branch cut down from a; = 0). The integral around 
the semicircle, 12, vanishes as the radius tends to infinity, 
while the integrals along 15/16 and 17/110 cancel exactly. 

The residue associated with the simple pole at = 1/2 
is cx exp(— it/2). Such constant amplitude oscillatory fre- 
quency contributions are typical of frequencies where there 
are zero-step solutions. Identical behaviour is found for the 
shell model, as discussed in Section |1] 

Working out the contributions from 13, 14, 18 and 19 
is more difficult. We will focus on 18/19, as the continuous 
spectrum contribution arises from the branch point at tj = 
X. We restrict our attention to the case when x is small (for 
illustrative purposes) . Making this assumption we note that 
ui and UJ — X are also small on 18 and 19. We will use Taylor 
expansions of the integrand using the three small parameters 
X, UJ, and uj — x . 

Keeping only terms up to second order, we find that the 
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contribution to equation 1B40I I from the integrals along 18 
and 19, lui, is 



ing the nature of solutions for large \uj\. We start by rewriting 
equation (I22I I as 



[ e-''^\Wo + Wiu;)\ai{x~uj)ln{x~uj) 

J 18.19 



I. = I e 

1 18.19 

+02(2; — uj)'^ ln{x — uj) + a3LLi{x — uj) ln{x — uj) 



+a4Uj ln(— cij) + a5Lo{x — lu) ln(— o;) 



(B41) 



where ai - as are constants, and we have expanded the 
inverse Wronskian in terms of the small parameter lu: 



Wiuj)'^ = Wo + WiLo + 0{lj^), 



(B42) 



where Wo and Wi are also constants. We now choose the 
phase of the logarithms such that 



ln(- 



ln(x' — uj) = 



ln\u;\ 
In la; 



2ni 



On 18 
On 19 



ln|a;-a;| On 18 

ln\x-Lu\+2TTi On 19. 



(B43) 



(B44) 



The logarithm terms cancel when we integrate over 18 and 
19, so the only non-zero contribution will come from the 27ri 
piece on 19. At first order, the contribution to is 



lull ~ l3i 



1 — ixt — e 



t2 



(B45) 



where /3i is a constant. The second order contribution to 1^^ 
is 



Iui2 = 



^ [71 + J^xt + ^-iX^t^ + (5l + 52Xt + S'iX^t^ 



(B46) 

where 71 — 73 and 5\ — Sj, are constants. If we include higher 
order terms we find that that the most persistent contribu- 
tion to U{x, t) arising from the continuous spectrum decays 
as 1/t. Higher order terms merely introduce contributions 
that decay as successively higher order powers of t. It is also 
clear that the continuous spectrum contributes a term to 
U{x,t) with a frequency that depends on position x. 

The analysis presented in this section is for a specific 
choice of initial data and for small x. We must therefore ask 
whether the key results would change if we had made differ- 
ent choices. It is certainly by no means clear that the contin- 
uous spectrum should decay for all smooth initial data. Al- 
though we have not observed anything other than decay, we 
cannot rule out the possibility of constant amplitude oscilla- 
tion or even growth. In contrast, the existence of a frequency 
component that depends on position appears inevitable. 



APPENDIX C: VANISHING OF SEMICIRCLE 
INTEGRAL IN DERIVATION OF INSTABILITY 
CRITERION 

In this Appendix we will demonstrate that we can neglect 
the integral over the semicircle in equation (I49II by consider- 



\/lu -^^n\/lu -n'uj ^0. (ci) 

We expand U in terms of the small parameter m/ui 



U = Uo + -Ui. 

UJ 



(C2) 



Equation ICH becomes 



VlUo ^ — V^f/i mVlUo + —Q.VlUi 

UJ UJ I UJ 

-n'Uo - —n'Ui] = 0. (C3) 

At leading order equation IIC3ll becomes 



V^[7o = 0. 



(C4) 



At first order equation ljC3|l becomes 



VlUi+n'Uo^Q. (C5) 
The general solution to equation ljC4^ is 



Uo = 



UoR = ci(f:pi)' 



x > y 



UoL = C2 f±i +C3 f-i - x<y, 



(C6) 



where for continuity of the function at 2; = y we require 



Cl = C2 



y-1 



+ C3. 



(C7) 



It is clear from equation IIC6t that Uo does not depend on 
UJ. Using the method of Green's functions to solve equation 
<U5ll we must first solve 



V^Go = S{x - x') 



(C8) 



The solution to this equation is given in Appendix^ equa- 
tions EUl and lEHl. We find Go (a; , x') to be independent of 
UJ. We now have 



UiR = - / Go{x,x')n' {x')Uor{x' ,y)dx' x>y 



UiL = - [ Go{x,x')Q'{x')UoL{x',y)dx' x<y. 

(C9) 

It is clear from equation ljC9^ that Ui must also be indepen- 
dent of UJ. Let us consider the implications of this for the 
Wronskian W . We get 

W = UoLU[,n ~ UorU;,!^ + -\UiLU;,n ~ UorU[l 

UJ L 

+UoLU[n ~ UirUL]+ O [^^y . (CIO) 
This means that we can write W = {1 — x^)W as 



W = fix,y) -I 9{x,y), 

UJ 



(Cll) 
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and we obtain 



W' ^ mg{x,y)/Lu^ ^ mg{x,y)f{x,y) 

W f{x,y) + mg{x,y)/uj ^ ui^ 

Integrating over frequency, 



W^^^^_ mg{x,y)fix,y) _ 
W 



In the limit as cj ^ oo the integral clearly tends to zero. We 
can therefore neglect the integral over the infinite semicircle 
in equation 1491 . 



